Alpha-particle condensation in ie O via a full four-body OCM calculation 
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In order to explore the 4a-particle condensate state in 16 0, we solve a full four-body equation of 
motion based on the 4a OCM (Orthogonality Condition Model) in a large 4a model space spanned 
by Gaussian basis functions. A full spectrum up to the 0^" state is reproduced consistently with 
the lowest six + states of the experimental spectrum. The Og state is obtained at about 2 MeV 
above the 4a breakup threshold and has a dilute density structure, with a radius of about 5 fm. The 
state has an appreciably large a condensate fraction of 61 %, and a large component of a+ 12 C(0^) 
configuration, both features being reliable evidence for this state to be of 4a condensate nature. 



It is well established that a-clustering plays a very 
important role for the structure of lighter nuclei 0, [4]. 
The importance of a-cluster formation also has been dis- 
cussed in infinite nuclear matter, where a-particle type 
condensation is expected at low density Q , quite in anal- 
ogy to the recently realised Bosc-Einstein condensation 
of bosonic atoms in magneto-optical traps 0]. On the 
other hand, for trapped fermions, quartet condensation 
also is an emerging subject, discussed, so far, only the- 
oretically [5j. In nuclei the bosonic constituents always 
are only very few in number, nevertheless possibly giv- 
ing rise to clear condensation characteristics, as is well 
known from nuclear pairing [6j. Concerning a-particle 
condensation, only the Hoyle state, i.e. the 0^~ state in 
12 C has clearly been established, so far. Several pap ers 
of the past 0, H, HI EU an d also more recently 11. 12, 11] 
have by now established beyond any doubt that the Hoyle 
state, only having about one third of saturation density, 
can be described, to good approximation, as a product 
state of three a-particles, condensed, with their c.o.m. 
motion, into the lowest mean field O^-orbit [14I. Il5j|. 

The establishment of this novel aspect of the Hoyle 
state naturally leads us to the speculation about 4a- 
particle condensation in 16 O, which is the focus of this 
work. The + spectrum of 16 has, in the past, very well 
been reproduced up to about 13 MeV excitation energy, 
including the ground state, with a semi-microscopic clus- 
ter model, i.e. the a+ 12 C OCM (Orthogonality Condi- 
tion Model) [l6| ■ In particular, this model calculation, as 
well as that of an a+ 12 C Generator-Coordinate-Method 
one [13], demonstrates that the 0^" state at 6.05 MeV and 
the 0^" state at 12.05 MeV have a+ 12 C structures [3 
where the a-particle orbits around the 12 C(0^~)-core in 
an S-wavc and around the 12 C(2^~)-core in a D-w&ve, 
respectively. Consistent results were later obtained by 
the 4a OCM calculation within the harmonic oscilla- 
tor basis [n|. However, the model space adopted in 
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is not sufficient to account simultane- 
ously for the a+ 12 C and the 4a gas-like configurations. 
On the other hand, the 4a-particle condensate state was 
first investigated in Ref. [ll( and its existence was pre- 
dicted around the 4a threshold with a new type of micro- 
scopic wave function of a-particle condensate character. 
While that wave function can well describe the dilute 
a cluster states as well as shell-model-like states, other 
structures such as a+ 12 C clustering are smeared out and 
only incorporated in an average way. Since there exists 
no calculation, so far, which reproduces both the 4a gas 
and a+ 12 C cluster structures simultaneously, it is crucial 
to perform an extended calculation for the simultaneous 
reproduction of both kinds of structures, which will give 
a decisive benchmark for the existence of the 4a-particle 
condensate state from a theoretical point of view. 

The purpose of this Letter is to explore the 4a con- 
densate state by solving a full OCM four-body equation 
of motion without any assumption with respect to the 
structure of the 4a system. Here we take the 4a OCM 
with Gaussian basis functions, the model space of which 
is large enough to cover the 4a gas, the a+ 12 C cluster, 
as well as the shell-model configurations. The OCM is 
extensively described in Ref. [20| . Many successful appli- 
cations of OCM are reported in Ref. @|. The 4a OCM 
Hamiltonian is given as follows: 



i i<j 
4 

+V£\i,j)] + ]T V 3a (i,j,k) + V ia (l,2,3A), (1) 

i< j<k 

where T h V^\i,j), v£?(i,j), V 3a (i,j,k) and 
V4 a (l, 2, 3, 4) stand for the operators of kinetic en- 
ergy for the i-th a particle, two-body, Coulomb, 
three-body and four-body forces between a particles, 
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respectively. The center-of-mass kinetic energy T cm is 

(p) 

subtracted from the Hamiltonian. V 2q is the Pauli 

exclusion operator [2l| . by which the Pauli forbidden 
states between two a-particles in 05, 0D and IS* states 
are eliminated, so that the ground state with the 
shell-model-like configuration can be described correctly. 
The effective a-a interaction is constructed by the 
folding procedure from two kinds of effective two-nucleon 
forces. One is the Modified Hasegawa-Nagata (MHN) 
force [12] and the other is the Schmidt- Wildermuth 
(SW) force 0, see Refs. [U] and 
respectively. We should note that 



24( for applications, 
the folded a-a 
potentials reproduce the a-a scattering phase shifts and 
energies of the 8 Be ground state and of the Hoyle state. 
The three-body force is phenomenologically introduced 
so as to fit the ground state energy of the 12 C. The 
same force parameter set as used in Refs. [HI or [24| 
is adopted in the present calculation. In addition, 
the phenomenological four-body force is adjusted to 
the ground state energy of 16 O. The three-body and 
four-body forces are short-range, and, hence, they only 
act in compact configurations. 

Utilizing the Gaussian expansion method [251 ] for the 
choice of variational basis functions, the total wave func- 
tion '3/ of the 4a system is expanded in terms of Gaussian 
basis functions as follows: 

*(0+)=5>^)$ c H, (2) 

$ c (i/) =S^pi- L (r 1 ,yi)ip h (r 2 ,y 2 )]i- L2 ipi 3 (ra,^3) , (3) 

where T\, r 2 and r% are the Jacobi coordinates describ- 
ing internal motions of the 4a system. <S stands for the 
symmetrization operator acting on all a particles obey- 
ing Bose statistics, v denotes the set of size parame- 
ters v\ , v 2 and 1/3 of the normalized Gaussian function, 
ifl(r, Vi) = Ni tUi r l exp (— Vir 2 )Yi m {r ), and c the set of rel- 
ative orbital angular momentum channels \\U , It\u^-, 1?^j 
depending on either of the coordinate type of K or H [25j | . 
where l\, l 2 and 1% are the orbital angular momenta 
with respect to the corresponding Jacobi coordinates. 
The coefficients A™ (is) are determined according to the 
Rayleigh-Ritz variational principle. 

Figure [T] shows the energy spectrum with J 7 ' = + , 
which is obtained by diagonalizing the Hamiltonian, 
Eq. |T]), in a model space as large as given by 5120 Gaus- 
sian basis functions, Eq. (J3]). It is confirmed that all 
levels are well converged. With the above mentioned ef- 
fective a-a forces, we can reproduce the full spectrum 
of + states, and tentatively make a one-to-one corre- 
spondence of those states with the six lowest + states 
of the experimental spectrum. In view of the complex- 
ity of the situation, the agreement is considered to be 
very satisfactory. We show in TABLE Q] the calculated 
rms radii and monopole matrix elements to the ground 
state, together with the corresponding experimental val- 
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FIG. 1: Comparison of energy spectra between experiment 
and the present calculation. Two kinds of effective two- 
body nucleon-nucleon forces MHN and SW are adopted (see 
text). Dotted and dash-dotted lines denote the a+ 12 C and 
4a thresholds, respectively. Experimental data are taken from 
Ref. [261 ] . and from Ref. [271 ] for the state. The assignments 
with experiment are tentative, see, however, detailed discus- 
sion in the text. 



TABLE I: The rms radii R and monopole transition matrix 
elements to the ground state M(E0) in units of fm and fm 2 , 
respectively. -R exp . and A/(E0) CX p. are the corresponding ex- 
perimental data. 
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ues. The M(E0) values for the 0J,0J, 
are consistent with the corresponding experimental val- 
ues. As mentioned above, the structures of the 0j and 
states are well established as having the a± 12 C(0 ] f ) 
and a± 12 C(2^") cluster structures, respectively. These 
structures of the 0^" and O3 states are confirmed in the 
present calculation. We also mention that the ground 
state is described as having a shell-model configuration 
within the present framework, the calculated rms value 
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agreeing with the observed one (2.71 fm). 

On the contrary, the structures of the observed O4 , 0^" 
and 0g states in Fig.[T]have, in the past, not clearly been 
understood, since they have never been discussed with 



K 



Al- 



the previous cluster model calculations 
though Ref. [ll|] predicts the 4a condensate state around 
the 4a threshold, it is not clear which of those states cor- 
responds to the condensate state. For example, the Ot 



state has been considered as one of the candidates 27 1, 
and also the state with large M(E0) value [H, 29. j, 
since the strong monopole transition implies a developed 
cluster structure [3(J. This ambiguous situation stems 
from the fact that the a-condensate-type wave function 
given in Ref. ll| can only reproduce the ground state 



and the 4a condensate state, but not the a+ 12 C configu- 
rations which make up a large part of the 16 0-spectrum 
up to the 0g state. 

As shown in Fig. [IJ the present calculation succeeded, 
for the first time, to reproduce the O4 , 0§ and 0g states, 
together with the Of , 0^" and 0^" states. This puts us in 
a favorable position to discuss the 4a condensate state, 
expected to exist around the 4a threshold. 

In Table HI the largest rms value of about 5 fm is found 
for the 0$ state. Compared with the relatively smaller 
rms radii of the O4 and 0^" states, this large size sug- 
gests that the 0g~ state may be composed of a weakly 
interacting gas of a particles [3l| of the condensate type. 




FIG. 2: 
with L 
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(Colors online) The radial parts of single-a orbits 
= belonging to the largest occupation number, for 



the ground and states with MHN force. 

While a large size is generally necessary for forming 
an a condensate, the best way for its identification is to 
investigate the single-a orbit and its occupation probabil- 
ity, which can be obtained by diagonalizing the one-body 
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(a) density matrix [14 



in the following way: 



J dr'p(r,r')^(r') = ^ L L (r), 



(4) 



where <fi L (r) is the single-a natural orbit with orbital an- 
gular momentum L, and fi L is its occupation probability. 
The one-body density matrix elements p(r, r') have been 



calculated as in 
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As a result of the calculation of 



the L = large occupation probability of 61% of 

the lowest OS-orbit is found for the 0g state, whereas the 
other five + states all have appreciably smaller values, 
at most 25% (0^). The corresponding single-a S orbit 
is shown in Fig. [2] It has a strong spatially extended 
behaviour without any node (OS*). This indicates that 
a particles are condensed into the very dilute OS single- 
a orbit, see also Ref. (33|. Thus, the 0^ state clearly 
has 4a condensate character. We should note that the 
orbit is very similar to the single-a orbit of the Hoyle 



state 
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/ery 

15J . We also show in Fig.[2]the single-a orbit for 
the ground state. It has maximum amplitude at around 
3 fm and oscillations in the interior with two nodal (2S) 
behaviour, due to the Pauli principle and reflecting the 
shell-model configuration. 
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FIG. 3: (Colors online) ry(r) defined by Eq. © for the 
state with the MHN force. 

In order to further analyze the obtained wave func- 
tions, we calculate an overlap amplitude, which is defined 
as follows: 

y(r) = ([^^^(r')<f>L( 12 C)]J*(0+)). (5) 

Here, $i( 12 C) is the wave function of 12 C, given by the 
3a OCM calculation [Hj], and r is the relative distance 
between the center-of-mass of 12 C and the a particle. 
From this quantity we can see how large is the component 
in a certain a+ 12 C channel which is contained in our 
wave function (2) for 0g . The amplitudes for the Og^ 
state are shown in Fig. [3] It only has a large amplitude 
in the a+ 12 C(0j) channel, whereas the amplitudes in 
other channels are much suppressed. The amplitude in 
the Hoyle-state channel has no oscillations and a long tail 
stretches out to ~ 20 fm. This behaviour is very similar 
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to that of the single-a orbit of the state discussed 
above. 

The a decay width constitutes a very important in- 
formation to identify the 0^ state from the experimental 
point of view. It can be estimated, based on the i?-matrix 
theory, with the overlap amplitude Eq. (O 34|. We find 



discussions concerning the 16 O spectrum with A. Richter 
are very much appreciated. 



that the total a decay width of the 0g" state is as small as 
50 keV (experimental value: 166 keV). This means that 
the state can be observed as a quasi-stable state. Thus, 
the width, as well as the excitation energy, are consistent 
with the observed data. All the characteristics found 
from our OCM calculation, therefore, indicate that the 
6th + state shall be identified with the experimental 0^ 
state at 15.1 MeV and that it is the a-condensate state 
of 16 0. 

Finally we discuss the structures of the 0^ and 0g 
states. Our present calculations show that the O4" and 
0^ states mainly have a+ 12 C(0^) structure with higher 
nodal behaviour and a+ 12 C(l~) structure, respectively. 
Further details will be given in a forthcoming extended 
paper. The calculated width of the 0^ is ~ 150 keV, 
which is much larger than that found for the 0^" state 
~ 50 keV. Both are qualitatively consistent with the cor- 
responding experimental data, 600 keV and 185 keV, re- 
spectively. The reason why the width of the 0|" state 
is larger than that of the 0^ state, though the 0|" state 
has lower excitation energy, is due to the fact that the 
former has a much larger component of the a+ 12 C(0^) 
decay channel, reflecting the characteristic structure of 
the Ot state. The 4a condensate state, thus, should not 



be assigned to the 0J or 0g state [35j but very likely to 
the 0^ state. 

In conclusion, the present 4a OCM calculation, for the 
first time, succeeded in describing the structure of the 
full observed + spectrum up to the 0<i~ state in 16 O. 
The + spectrum of 16 up to about 15 MeV is now es- 
sentially understood, including the 4a condensate state. 
This is remarkable improvement concerning our knowl- 
edge of the structure of 16 O. We found that the 0^ state 
above the 4a threshold has a very large rms radius of 
about 5 fm and a rather strong concentration of a parti- 
cles 61% on a spatially extended single-a OS* orbit. The 
wave function was shown to have a large a+ 12 C ampli- 
tude only for 12 C (Hoyle state). These results are strong 
evidence of the newly found 0^ state as being the 4a 
condensate state, i.e. the analog to the Hoyle state in 

C. Further experimental information is very much re- 
quested to confirm the existence of this novel state. Also 
independent theoretical calculations are strongly needed 
for confirmation of our results. 
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